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The finite element technique is used widely for researchers and 
manufacturers to design and simulate electrical systems in general and 
electrical machines such as shunt reactors (SRs) and transformers in 
particular. Many papers have recently applied several methods to analyze 
magnetic fields, copper losses and joule power losses in the shunt reactors 


(SRs). In this research, the finite element technique with coupling to global 

quantities is proposed to investigate the voltage and current distributions in 
Keywords: the windings, and compute the distribution of magnetic field in the air gap 
and along the air core of the SR, as well as copper and core losses. The 
developed method is directly applied to the practical SR of 91 MVAr and a 
Copp er loss rated voltage of 500 kV. The finite element method (FEM)-simulated results 
Finite element method are validated with experimental results to ensure accuracy and reliability. 
Iron loss This facilitate designing the reactor. 
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1. INTRODUCTION 

Nowadays, the shunt reactor (SR) is widely applied to stabilize the electrical systems due to 
parasitic capacitance appearing in the transmission line. It is also applied to bound the short-circuit current 
and overvoltage in the electrical system. In addition, it can be also used to imbibe not only the reactive power 
created by a conductive capacitance in power transmission lines, but also adjust the reactive power in the 
electrical system. Recently, several studies have presented the many different method to investigate and 
analyse electromagnetic parameters of the SR [1]-[10]. In reference [1], the finite element method (FEM) 
was presented to investigate the flux fringing distribution around the air gap of the magnetic circuit of the 
SR. In reference [2], the theory of Maxwell stress tensor was proposed to evaluate the electromagnetic force 
on the core blocks causing the vibration and noise during the working of SRs. In reference [3], the paper 
investigated the influence of distance air-gaps on the iron-core of SRs by using the analytical method. In this 
work, the winding inductance was also computed to define the relationship of the distance between the air 
gaps. In references [4], [5], a mathematical model was studied to consider the nonlinear dynamic case. In 
references [6], [7], the relation between the SR and power system in various transient situations was 
investigated. In reference [8], the effects of core block numbers on the inductance was studied via the 
analytical technique. In reference [9], a FEM was presented to compute the magnetic flux density distribution 
and losses in the magnetic circuit of the SR. In reference [10], [11], the papers provided a useful tool for the 
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design and optimization of ferromagnetic-core inductors, which are commonly used in a variety of electronic 
devices. The analytical model presented can help engineers to reduce the time and cost involved in designing 
these devices while ensuring that the design meets the required specifications. In reference [12], this paper 
provides valuable insights into the use of air gaps in power inductor design. By understanding the influence 
of air-gap arrangements on inductor performance, designers can create more efficient and reliable inductors 
for a variety of applications. 

It can be seen that the above studies have not investigated the voltage and current distributions on 
the windings, magnetic flux density on core blocks and in air gap between core blocks. In addition, the 
inductances of each winding phase has been not also considered in the previous studies. In this research, a 
FEM is developed to compute and design an oil-immersed SR of 91 MVAr, 50 Hz, 500 kV. The scenario of 
this paper is divided into two senarios. In the first approach, an analytic model is developed to define the 
desired output parameters of the proposed SR. In the second approach, a FEM is developed with the obtained 
parameters in the first step to compute and analyse the voltage and current distribution, magnetic flux density 
distribution, Joule power loss density and winding inductances. The simulated results of the FEM are verified 
with the experimental results to validate on the proposed method. 


2. ANALYTICAL APPROACH OF THE SR 

One of the most important things in designing the SR is the factor of the air gap volume. This factor 
depends on parameters of the SR, i.e. reactive power, magnetic flux distribution in the magnetic circuit, 
winding inductances, power frequency, and energy stored in the air gaps and the windings. For the SR, the 
reactive power (Q) is determined as (1) [1], [10], [13]. 


Q= Uef lef (1) 


Where Ues and Iep are respectively the effecitive voltage and current. It should be noted that the resistance of 
windings is small in comparison with inductances, therefore this value can be ignored. For that, the 
electromagnetic force (EMF) and effective current are determined as (2) and (3) [13]. 


2T 
Eef = Ues = (=) FN Ba Ape: (2) 
— Rgap _ (1)\ 3Bm-gap 
I= N (=). Uo-N (3) 


Where f is the frequency, N is the turn number, B,, is the maximum flux density, Agap is the area of the air 
gap, [gap is the length of the air gap, Ho is the air permeability and Rgap is the air gap reluctance. Based on 
(2) and (3), the air gap volume V,,,, can be presented via the air gap length lgap: 


l =n (4) 


gap © Tepe 
mo me 


Vjap = Agap: 


the air gap volume is considered as a constant if the quantity Q, magnetic flux and frequence are the constant. 
Therefore, the dimension of core (D,) is determined as (5) [14], [15]. 


= [+s (5) 


The inductance of winding is calculated via the reactive power and voltage as (6). 


u2 
— ef 
L=% (6) 


The size of the SR is influenced by the air gap length and area. Thus, the iron losses, copper losses 
and other factors depend on the different ratio value of the area and air gap length. Because this factor will 
decide to the cost and losses of the SR. The air gap permeance can be defined as (7) [16]. 


Aga 
Pray = (222) uo = P, (7) 


lgap 


The relation between the winding inductance and air gap permeance is: 
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the winding inductance can also be computed via the expression as (9) [10], [13]. 
2 
L= (=) (9) 


From (8) and (9), one gets (10). 


E 
N= |ia = A (10) 
tolga? [wf po GED 


3. FINITE ELEMENT APPROACH 
3.1. Canonical magnetodynamic model 

A model of a studied problem with a studied domain 0 = NQ, U Q% (where domain 0, is the 
conducting region and and Q¢ is the non-conducting region) is presented in Figure 1. The boundary of Q 
denotes 02 =I = fj, U Ik). The set of Maxwell’s equations and behavior laws are written in the frequency 
domain as [10]: 


rot H = J, (11) 
rot E=—jwB (12) 
div B =0 (13) 
B = uH (14) 
J=0E (15) 


where H: the magnetic field (A/m), B: magnetic flux density, E: electric field (V/m), J,: electric current 
density (A/m’), J: eddy current density (A/m?), u: relative permeability, and ø: electric conductivity (S/m). 


T=T,uUl. 


Figure 1. General magnetodynamic model [10] 


These fields defined in function spaces (H, (rot; 2) and H, (rot; 2) must satisfy Tonti’s diagram 
[10]. These function spaces contain the fields defined on I} and I, of the studied domain and BCs on I as 
(16) and (17) [10], [17]-{20]. 


nx H|p, = 0 (16) 

n+ Bir, =0 (17) 
Where n is the unit normal exterior to Q. The set of (11), (12), and (13) is solved with the BCs in (14) 
and (15). In addition, global quantities with the current (/;) and voltage (V;) imposed the inductor are shown in 


Figure 2. The source of electromotive force (EMF) is defined through the surface Iji» i.e. (18) and (19) 
[8], [21]. 
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PE 

GE dl = v; (18) 
E 

$èinJds=h (19) 


Figure 2. Magnetodynamic model with global quantities (current and voltage) [10] 


3.2. Magnetic vector potential weak formulations 

The magnetic vector potential weak formulation (A) (B = curl A and E = gð ,A — ograd v) is 
defined via the Ampere’s law (10 A) in Q [10], [22]-[24]. 

(curl A, curl w’) + (oð; A,w')a, + (ograd v, w')a, +< n X H,w' >p,-r, 

Q 

+< [n x H]r w >r = Uwa p Vw' E H+ (O) (20) 
Notations of (.,.)o and <., >r, are the volume and surface integrals, respectively, where the volume integral 
is defined in Q and the surface integral is defined on T. The function space H+ (N) contains the test function 
w’ and basis functions for A. At the discrete level, this function is determined by edge finite elements (FEs), 
which are a type of numerical discretization technique used to approximate the solution of partial differential 


equations. The boundary < n x H, w’ >r,-r, in (20) represents the natural BCs of (15), typically zero. From 
(20), by taking w’ = grad v’ as a test function in the weak formulation (15), one gets (21). 


(00,A, grad v')o, + (ograd v, grad v')a, +< [n X hi]p,w’ >r 
=< n: j, v >r, Yw' E H! (N) (21) 


Where I, is one of the surfaces I ;, which is fixed by a current or voltage as shown in Figure 2. 


3.3. Coupled to voltage quantity 

The voltage applied (V;) to the inductor can be expressed as a natural global constraint via an unit 
source electric scalar potential v,; associated with the applied voltage. The voltage applied to each massive 
inductor (Q,;) is expressed as a sum of the applied voltages between the electrodes i; and I, . This results 
in an equation expressed as the circuit relation associated with the inductor Qm, via the voltage V; [4], [25] 


vi = ViVivs i (22) 
By substituting (17) into (16), one gets (23). 


(curi A, curl w’) + (oð,A,w')o, + Xi Vi (ograd vsiw Jas +< nX H, w >rp-r 
a 
+< [nx H] w >n= U,w')o,Vw' E HEO) (23) 


3.4. Computation of Joule power losses via a post-processing 
Joule losses are then computed with [21], [22]: 


1p JI 
Pioss = zho zan (24) 


where J is the conjugate of J. By using the Poynting theorem associated to the surface integration of degrees 
of freedom of A and v located on the border, the volume integration can be defined as (25). 
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Poss = ŻRe ($ (nx H) dr) = Re (-“¢ (n x H, —n xX grad v) 4dr) (25) 
loss ~ 5 T ~ 2 Ir S g ; 


Where E and A are respectively the complex conjugates of electric field E and magnetic vector A. 


4. APPLICATION TEST 

The test problem is herein a three phase SR as pointed out in Figure 3, where Figure 3(a) illustrates 
the representation of the air gaps between core blocks, Figure 3(b) describes the winding structure and 
Figure 3(c) depicts the structure of a three phase SR with the winding and magnetic circuit. The main 
parameters given in Table 1. The 3D model of the SR is established as in Figure 4, where a full 3D model is 
given in Figure 4(a). Due to the symmetric nature of the structure of the SR, in order to reduce computational 
time, a half cross-sectional model is considered in Figure 4(b), and imput parameters for constructing the 
model are pointed out in Figure 4(c). 


Ww ol Lew Li Dy =! La We m 


(a) (b) (c) 


Figure 3. Geometry of a three phase SR (a) air gaps between core blocks, (b) winding structure, and 
(c) construction of a a three phase SR 
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Figure 4. 3D-model of a three phase SR (a) 3D model of SR, (b) a half cross-sectional model, and 
(c) input parameters 
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Table 1. Parameters of the SR 


Parameter Symbol Value Parameter Symbol Value 
Reactive power Q(MVAr) 91 Deep Yoke Dy (mm) 776 
Rated current I(A) 105 Total air gap length along the core l (mm) 400 
Rated voltage U (kV) 500 Number of turns N (vòng) 2244 
Inductance Xı (Q) 2712 Distance between the core and windings bew (mm) 135 
Core diameter D, (mm) 666 Yoke width Wwy (mm) 236 
Core height H. (mm) 1978 Winding height Hy (mm) 1559 


The distribution of phase voltage in the windings is shown in Figure 5. It shows that the effective 
voltage applied to windings (phase A, phase B, and phase C) is the rated phase voltage with the value of 
288,675 kV, which is close to the measured result of 288,7567 kV. This means that the significant error on 
the voltage simulation is less than 1%. 

In the same way, the phase current distribution in the winding is presented in Figure 6. The 
comparison of current and voltage values between the simulated and measured results is given in Table 2, 
with errors being lower than 1% for both cases. The magnetic flux density distribution in the magnetic circuit 
is presented in Figure 7. It can be shown that the value is the biggest when the current of phase B is the 
maximum. The magnetic flux density distribution along the line of C1-C2 (core block) is presented in 
Figure 8. It shows that the magnetic flux mostly focuses on the surface of the core blocks, thus the magnetic 
flux density near edges of the core blocks is bigger than the one inside core blocks. In order to increase the 
leakage flux in the vicinity of the core blocks, the reluctance of magnetic circuit and number of air gap core 
blocks need to increase also. 
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Figure 6. Current distribution in winding 
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Table 2. Comparison of simulated and measured voltages 


Parameters Rated values Simulated results Measured results Error (%) 
Current of phase A (A) 105 105.67 106.44 0.73 
Current of phase B (A) 105 105.71 106.86 1.08 
Current of phase C (A) 105 105.65 106.34 0.65 
Voltage of phase A (V) 105 288.62 288.76 0.35 
Voltage of phase B (V) 105 288.56 288.71 0.5 
Voltage of phase C (V) 105 288.62 288.93 0.12 


B Liew ted 


Figure 7. Distribution of magnetic flux density in the magnetic circuit 
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Figure 8. Distribution of magnetic flux density along the line of C1-C2 (core block) 


The magnetic flux density distribution along the air gap (G1-G2) between core blocks is indicated in 
Figure 9. The results show that the magnetic flux density along the inner winding surface is very small and 
uniform due to the large number of air gaps. Its mean value is equal to 0.1797 T. The values of self and 
mutual inductances between phases (phase A, phase B, and phase C) are given in Table 3. The comparison of 
phase reactances between the rated and simulated results is given in Table 4. It can be seen that the error is 
lower than 1% for phase B and equal to 1% for phase A and phase C. In the same way, Table 5 shows the 
comparison of simulated and measured reactances for three phases. The errors on them are smaller than 2%. 
This is a good agreement with the developed method. In particular, it is satisfied the IEC 60076-6 of the SR. 
The copper and iron losses of the SR are pointed out in Figures 10 and 11, respectively. The results of iron 
and copper losses in the SR between the simualtion and measurement are given in Table 6. The error is lower 
than 4% for the iron loss and 2% for the copper loss. These errors are completely acceptable and realiable. 


Table 3. The values of self and mutual inductances (H) 
Inductance (H) Phase A__ Phase B Phase C 


Phase A 8.7257 -0.00565 -0.00129 
Phase B -0.00565 8.7199 -0.00566 
Phase C -0.00129 -0.00566 8.7268 
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Table 4. Comparison of simulated and rated reactances 


Parameters Rated value (ohm) Simulated result (ohm) Error (%) 
Reactance (phase A) 2,712 2,739.08 1.00 
Reactance (phase B) 2,712 2,735.88 0.88 
Reactance (phase C) 2,712 2,739.42 1.00 


Table 5. Comparison of simulated and measured reactances 


Parameters Simulated result (ohm) Measured result (ohm) Eror (%) 
Reactance (phase A) 2,739.08 2,698.7 1.47 
Reactance (phase B) 2,735.88 2,694.4 1.52 
Reactance (phase C) 2,739.42 2,689.3 1.83 

0.80 
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Figure 9. Magnetic flux density along the air gap (G1-G2) between core blocks 
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Figure 10. Copper loss in the winding of the SR 
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Figure 11. Iron loss in the magnetic circuit of the SR 
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Table 6. Comparison of simulated and measured copper and iron losses 
Parameters Simulated result (kW) Measured result (kW) Error (%) 
Iron loss 68.519 70.817 3.24 
Coppeer loss 152.377 155.321 1.90 


5. CONCLUSION 

In this research, the proposed SR was designed using an analytical approach to determine its main 
parameters. Afterward, a FEM simulation was conducted to simulate the prototype of the SR. The results of 
the SR obtained from the FEM were then compared to the that from the measured method to confirm their 
accuracy. The obtained results have also suggested that the analytical approach and FEM simulation are 
accurate methods for calculating various parameters of the SR. These parameters include voltage, current, 
inductance, reactance, copper loss, and iron loss. Additionally, these methods can also evaluate material 
costs, losses and other critical parameters. Moreover, the use of FEM simulation allows for a faster design 
and computing process and to optimize the prototypes needed to achieve the expected results. This approach 
can potentially save time and cost in the design and manufacturing of shunt reactors. 
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